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ABSTRACT 


This paper derives (§§II-TV) and summarizes (§VI) a new "plug-in- 
and-grind" formalism for calculating the gravitational waves emitted by 
any system with weak internal gravitational fields. If the internal 
fields have negligible influence on the system’s motions, then the new 
formalism reduces to standard "linearized theory". Whether or not grav- 
ity affects the motions, if the motions are slow and internal stresses 
are weak, then the new formalism reduces to the standard "quadrupole- 
moment formalism" (§V). In the general case the new formalism expresses 
the radiation in terms of a retarded Green's function for slightly curved 
spacetime, and then breaks the Green 1 s-f unction integral into five 
easily understood pieces: direct radiation , produced directly by the 

motions of the source; whurap radiation , produced by the "gravitational 
stresses" of the source; transition radiation , produced by a time-changing 
time delay ("Shapiro effect”) in the propagation of the nonradiative, 1/r 
field of the source; focussing radiation , produced when one portion of the 
source focusses, in a time-dependent way, the nonradiative field of 
another portion of the source, and tail radiation , produced by "backscatter" 
of the nonradiative field in regions of focussing. 


I. INTRODUCTION . 7 

a) Introduction to This Series of Papers 

Thanks to the pioneering work of. Joseph Weber (1960,1969), "gravi- 
tational-wave astronomy" may be a reality by 1980. Although Weber’s "events" 
may turn out to be non-gravitational in origin, second-generation detectors 
of the Weber "resonant-bar" type, with amplitude sensitivities roughly 100-fold bet- 
ter than today's bars, are now under construction (Braginsky 197b, Fairbank and 
Hamilton, as described in Boughn et al . 197k); and third- genera tion detectors are 
being discussed. The third generation should be able to detect and study the 
gravitational-wave bursts generated several times per year by supernovae in the 
Virgo cluster of galaxies. Detectors with other designs may succeed in 
detecting waves from pulsars [see, e.g., Braginsky and Nazarenko (1971)] 
and from near-encounters of stars in dense star clusters [gravitational 
bremsstrahlung ; see. e.g., Zel’dovich and Polnarev (1974)3. 

And, of course, totally unexpected sources may be detected. [For reviews 
of the prospects for gravitational-wave astronomy see Misner (197k), Rees 
(197k), and Press and Thorne (1972).] 

In preparation for the era of gravitational- wave astronomy, our 
Caltech research group has embarked on a new project: We seek (1) to elu- 
cidate the realms of validity of the standard wave-generation formulas; 

(2) to devise new techniques for calculating gravitational-wave generation 
with new realms of validity; and (3) to calculate the waves generated by - 
particular models of astrophysical systems. Throughout this project we 
shall confine ourselves to general relativity theory. 

Most past calculations of gravitational-wave generation use one of 
three formalisms: (1) "linearized theory" or its quantum-theory analogue; 
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(2) the "quadrupole-moment formalisin''; (3) "first-order perturbations, 
of stationary, fully relativistic spacetimes." 

"Linearized theory" is the formalism obtained by linearizing gen- 
eral relativity about flat spacetime [see, e.g., Chapters 18 and 35 of 
Misner, Thome, and Wheeler (1973) — cited henceforth as "MTW"]. It is 
also the unique linear spin-two field theory of gravitation In flat space- 
time — and as such it has a simple quantum-theory formulation. [For ref- 
erences and overview see, in MU'/, §7.1, Box 7.1, and part 5 of Box 17.2] . 
Linearized theory is typically used to calculate wave generation when the 
source's self gravity has negligible influence on its motions (e.g., waves 
from spinning rods and from electromagnetic fields in a cavity) . In this 
paper we shall devise a new wave-generation formalism valid for any system 
with small but non-negligible self gravity; and in Paper III (Kovacs and 
Thorne 1975) we shall use that formalism to calculate the gravitational 
br eras str alllung produced when two stars fly past each other with large im- 
pact parameter, but with arbitrary relative .masses and velocities. 

The "quadrupole-moment formalism" (in which the wave amplitude is . 
proportional to the second time derivative of the. source's mass ’ quadrupole 
moment) dates back to Einstein (1918), and has been canonized by Landau 
and Lifshitz (1951). The derivations of this formalism which we find in 
the literature are valid only for systems with slow internal motions and 
weak (but non-negligible) internal gravitational fields [see, e.g., the 
past "^-Newtonian derivation by Chandrasekhar and Esposito (1970), the 
matched-asymptotic-expansion derivation by Burke (1971), and the de Donder 
gauge derivation by Landau ana Lifshitz (1951) as made more explicit in 
Chapter 36 of MTW] . However, a detailed analysis given in Paper II 
(Thorne 1975) shows that only the slow-motion assumption is needed: the 
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quadrupole-moment formalism is valid for any slow-motion system, regard- 
less of its internal field strengths. Paper II also extends that formal- 
ism to include the radiation produced by all of the source's other 
moments (both "mass" moments and "current" moments); and it derives 
formulas in terms of the moments for the near-zone fields, the radiation 
field, the radiation reaction, and the energy, momentum, and angular 
momentum carried off by the waves. In a forthcoming paper Thorne and 
Zytkow (1975) will use the extended formalism of Paper XI to calculate 
the "current-quadrupole" gravitational waves produced by torsional oscxl- 

lations of neutron stars. 

"First-order perturbations of stationary, fully relativistic space- 
times" is a technique that has been used extensively in recent years to 
analyze waves ! from "fast-motion" oscillations of black holes and neutron 
stars, and from particles moving In the Schwarzschild and Kerr gravita- 
tional fields. [For reviews, see Press (1974), Ruff ini (1973) and §36.5 
of MTW; see also the recent paper by Chung (1973).] It is not yet clear 
whether our project will delve into this technique. • 

b) Overview of This Paper 

In this paper we confine attention to systems with weak internal 
gravitational fields. Section II rewrites the exact Einstein field equa- 
tions in a non-covariant form ("de Bonder form") that is amenable to weak- 
field approximations. Section III gives a systematic account of approxi- 
mate, weak-f ield formalisms based on the exact de Donder form of the field 
equations — including the accuracy of the various formalisms and their re- 
lationships to each other. Section IVa applies the analysis of Section 
III to astrophysical systems, and concludes that, when analyzing their 
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structure and evolution, one. must typically calculate the stress-energy 
tensor and gravitational field with accuracies: 

| (error in 2 T ’~^ / 2 T °°I ^ 

| (error in /TP® | . £ £ 

£ 5 (typical value of inside source) 

n. (mass of source) / (size of source). 

Section IVa also concludes that the external gravitational field ^h^ 
must typically be calculated to accuracy . 

((error in 2 h^ V ) / 2**°° ^ ^ 

if one desires reasonable accuracy in the radiative part of that field. 

Section IVb presents a "post-linear" formalism for calculating a 
system’s structure and evolution (^T^ V and ^h^ V ) to the desired accuracy; 
and Section IVc derives a plug-in-and-grind formula for the higher- 
accuracy external field ( ? h ), which contains the radiation. Section V 

shows how the resulting formalism, when applied to slow-motion systems, 
reduces to the standard "quadrupole-moment formalism". ' 

We recommend that, before tackling the rest of this paper, the 
reader peruse Section VI. That section summarizes our post-linear for- 
malism and our formula for the external (radiation) field. 

The "guts" of this paper, in terms of complex calculations, reside 
in the Green ’ s-funct ion manipulations of Section IVc. Our particular way 
of. handling the Green's functions is motivated in Appendix A, and has been 
influenced by the following papers: De Witt and Brehme (1960) (exact Green’s 


4 



functions for scalar and vector wave equations in curved spacetime) ; 
Robaschik (1963) (exact Green's function for tensor wave equation in 
curved spacetime); John (1973a, b), Bird (1974) and especially Peters (1966) 
(Green’s functions in weakly curved spacetime). Although these papers had 
much influence on us, our specific manipulations are so different that we 
have found it impossible to trace the details of that influence in our writeup 

II. EXACT GENERAL RELATIVITY, REWRITTEN IN "WEAK-FIELD LANGUAGE" . 


We begin by writing the exact, nonlinear Einstein field equations in 
an arbitrary coordinate system in the form [§20.3 of MTW; §100 of Landau 
and Lifshitz (1962)] 


< -f ,0*6 ■ l^C-8)(T yv + t£ L ) 


where "L-L" means "Landau-Lifshitz” , and where 


, T yav3 „ • yv a3 av y3 . 
H L-L = 9 9 “ 9 9 


( 1 ) 


(2a) 


r 


c-g) “ ’" det llg uv li = -det||a yv || ; 


a3 - maw „ m- 1 / <*3 Ay aX 3y 
: L _ L = U6ir(-g)] {q >A q , v " 9 ,\9 

, ctX 3v yp , 3A av yp 
^ 8 s yv9 ,p9 ,X 8 8 yv9 ,p9 

,1 /0 aX 3y «3 Ay. , . 1 

+ 8 (2g g g g )( 2 S v p 8 aT ~ 8 ocj s vt^ 3 


(2b) 


. 1 a3 Av du 


+ 2 g 8 V,9 >P 9 ; 


* , vp aX 3 u 

X> + 8 Xy 8 9 I ,v3 ,p 

(2c) 

* v.r oa i 

vr 3 ,X 9 ,y ' 

(2) 


,yv 


The equations of motion for the material stress-energy tensor T follow 
directly from the field equations (1) and can be written in the equivalent 
forms 

r-jQ niRn ft 

(3) 


T a3 « 0 

. ;3 ’ 


[(-g)(T a3 + t“^ L )] >3 = o 
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[Here and throughout this series of papers we use the notation and sign 

oc 3 

conventions of MTW; in particular c = G = 1; and g are the com- 

ponents of the metric; commas and 3*s denote partial derivatives, Y a = 

3 Y = 3Y /3x Ct ; semicolons denote covariant derivatives with respect to the 

a 

metric g a g; and our signature is (-+++).] 

Now, and henceforth in this paper, we impose three restrictions on 
our analysis: (1) We confine attention to systems with "weak internal 
gravity" — i.e., to systems throughout which one can introduce nearly Lorentz 
coordinates. (2) We confine ourselves to "isolated systems" — i.e., to sys- 
tems that are surrounded by a region ("local wave zone") , much larger than 
a characteristic wavelength of the emitted waves, in which all waves are 
outgoing and in which external masses have negligible influence on the 
gravitational field. (3) We restrict our analysis to the interior of the 
source and its local wave zone, and throughout these regions we use nearly 
Lorentz, asymptotically flat coordinates, specialized to satisfy the deDonder 
gauge condition. 

Mathematically, these restrictions state that the "gravitational 
field" •• : 


e - 3 ^ + n yv O) 

has the properties 

jh :JV | << 1 everywhere, (5a) 

[h UV j ^ 1/r as r -* where r = (x 2 + y 2 + z 2 )"^ 2 , (5b) 

"h^ V is devoid of incoming waves at r -*■ 00 , (5c) 

h^ V =0 .(deDonder condition). (5d) 

» V 


With these restrictions, the exact Einstein field equations (1) take on the 


6 



form 


□ „ - -16ir(-g) 1/2 (T PV + t^ L )- (-g)" 1/2 t? Ja ft h'" 3 „ • (6) 


Here [3 , is the wave operator for scalar fields in the curved spacetime' 


described by the metric g 


ct3 


□ s = (~s)" 1/2 ((-g) 1/2 g a6 a ) . 


(7) 


[Appendix A, which is best read after one has finished reading the rest of 
the paper, explains why we write the field equations in terms of d 

s 

□ o* B 

f = n 

(the flat-space wave operator) or LH (the curved-space wave operator for 
tensor fields).] 

Equations (2b)- (7) are the exact, nonlinear equations of general 
relativity for any isolated, weak^fie Id system — but they are written in a 
very special coordinate system rather than in generally covariant form. 

Because |? JV | « 1, we can express each quantity in our formalism, 
except T^ V , as a power series in h^ V . When writing down such a power 
series, it is convenient to raise and lower indices of. with the 

Minkowski metric r}^ - T) = diag (-1,1,1,!) : 


h v = n 

a 'ay 


h a£5 = WW 1 , 


h = h 


a 


a 


etc. 


(8a) 


It is also convenient to define a "trace-reversed" gravitational field 


by 


h yv E h 1 ^ - i h rT , 


(8b) 


and to raise and lower its indices, like those of h : , with the Minkowski 

metric. Note that equation (8b) implies 
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. _ . a 7- -ruv . pv 1 , mv 
h = h = -h , h- ~ h H - -7- h ry 

a ■ 2 


(8c) 


To derive the explicit power series expansions for g^ V , g^ v » 
etc. , one can proceed as follows. Equation (4) is the desired expansion 
for . It contains only two terms: 


q .MV = n MV _ 


(9a) 


The expansion for the metric determinant (-g) is obtained by inserting ex- 
pression (9a) into the second of equations (2b): 


(~g) = -det |jg 1JV || = -det||n 1JV “ h^H 


1 - h + \ [(h) 2 - h 01 ^] +0[(h) 3 ] 


(9 b) 


The contravariant components of the metric are' then obtained by inserting 
(9a, b) into the first of equations (2b) 


mv 


c-g) 9 


- n pv - h h J ' J + | n uv [ (h) 2 + 2^ B h„ D ] +o][(h) 3 ] 


ag- 


(9c) 


and the covariant components are obtained as the matrix inverse of these 
contravariant components: 


g = n +11, -|hn + h h® - ~ h h + -~ n ,(h 2 - 0 )+ or (h) 3 l 

6 MV mv yv 2 yv \ia v 2 mv 8 yv v ag^ UlLvn; J 


(9d) 

The connection coefficients , which appear in the usual expression 

. = 0 for the. equations of motion, are obtained by inserting expansions 

5 V 
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(9c, d) into the standard formula 


f y _1 UY/ , _ v 

a3“2 g ^ s ya ,3 g y3,a 8 a3,Y ; 


“ "Xx YfXg h ,a- V h ’ V)+0[(h)2] (9e) 


Similarly, the scalar-wave operator Q is obtained by inserting expan- 

s 

sions (9b, c) into equation (7), and using the deDonder gauge condition 
(5d) to simplify: 


□ - n a3 a Jo - “ 4 h n ap )3 J c + 0 [(h)^] 


^3- 1 TT a3. 


a 3 


a 3 


(9f ) 


and the components of the Landau-Lif shitz pseudotensor are obtained by in- 
serting expansions (9a,b,c,d) and the deDonder condition (5d) into equa- 
tion (2c) : 

* ■ a .-.- 1 ii <V.. * v* l6u .<, 


- (n“\ 5*” 5» . + n s \i , E™ „ V , > 

' ,p ,x . fp fA 


( 9 g) 


* | < 2 „“V - »“V“X2 V „ - .„V S"- • 


IJV ^ v -tjlV 

Henceforth in this paper we shall regard q , (~g) > g , g^, 

r y _ T~1 . and t?^_ as shorthand notation for the infinite power series 
aS L-L . ~ 

expansions, whose first few terms are shown in equations (9), Given these 
expansions, the full content of general relativity is embodied in the equa- 
tions of motion for the matter 
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If the system is several lumps, with masses m and sizes ft , separated 
by distances b » ft (e.g., a binary star system or two stars flying past 
each other) , then ' 

£ ^ m/b if one is interested only in the relative motions of the lumps 

£ ^ m/ft if one is also interested in the internal structure and 
dynamics of the lumps. 

We shall characterize every weak-field approximation formalism by 

two integers and . These "order indices" tell us the magnitude 

2 ' ' 

of the errors made by the .'formalism : 

^Note that all of the | T ; “ ^ | are T^, and consequently all of the 
|h^ V | are £ h . This fact dictates the form of equations (12), 


n„ 


(errors, in T 1JV ) / T°°| ^ e T 
(errors in h^ V ) / h^j r 


(12a) 

(12b) 


For example, a formalism of order (n^n^) = (l,l) makes fractional errors 
of order £ in both the stress-energy tensor and the gravitational field, 

while a formalism of order (2,1) makes fractional errors z~ in and 

■rMV '■ 

£ xn h 

Errors in when fed into the equations of motion (l'Oa), produce 

errors in and similarly, errors in T^ V , when fed into the field equa- 

tions. (10b), produce errors in h^ V . This feeding process places con- 
straints on the order indices (n_,n, ) of any self-consistent approximation 

T h 

formalism. The constraints are revealed explicitly by an order-of-magnitude 
analysis of equations (10a, b): 
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Consider a weak-field system with characteristic field strength 
e and characteristic length-time scale & . Below each term of equations 
(10a, b) , write the order of magnitude of that term: 






„ r d T ctv _ p v T ya 
av av 


(13a) 


(T 00 /Jt.) (e/«(T 00 ) (e/S,)(t 00 ) 


=-16*(-g) 1/2 I^ V -16w(-s) 1/2 t^ L -C-g)' 1/2 ^ a >8 ^ 8 >a • (13b) 


(e/S 2 ) 


,00 


2 2 

(e z /JT) 


7 7 

(e‘ l /tn 


^ — »JV 

Equation (13a) shows that fractional errors £ in h’ produce frac- 


tional errors £ 


V - in T^; i.e., i.e., 


n T i "h + 1 


(14a) 


Equation (13b) — together with the order-of-magnitude field equation 

00 2 n\) 

T *V/ e/5/ — shows that fractional errors e in T“ produce fractional 

errors £ ^ in ; i.e., >_ £ i.e., 


< n T . (14b) 

Equations (14a, b) can be restated as the following constraints on the order 
indices of any self-consistent approximation formalism: 

n h = n T - 1 or n,, = n. f , . (15) 

In other words, the order (n^,n^) of any approximation formalism must be 
either. (n,n~l) or (n.n) for some integer n . 

Suppose that a specific system has been analyzed using an approxima- 
tion formalism of order (n,n-l). Denote by (x ) and ( x ) c ^ e 
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explicit expressions obtained in that analysis for the. system's stress- 
energy tensor and gravitational field* From these expressions it is 
straightforward to generate an "improved" gravitational field ^h (x ) 
with fractional errors £ n . The key to doing this is the structure of 
the field equations (10b): In these field equations, fractional errors 

of order ^ in produce fractional errors of order £ U in both 

G and the expression 

S' - 

- ifiuc-g) 172 d wv + O - (-g)- 1/2 h^ >B • 


Hence, ^h^ V (t,x) satisfies the differential equation 

, t .D P" - , n . . C-167TC-S) 1/2 (X^ + 1^\) - (-g)- 1/2 g^ 6 J . 

(n,n-l) s n (n,n~l) u ° L-L »p ,a 

(16) 


Here the prefix (n,n-l) means that a quantity is to be calculated, with 
fractional error £ n , using and • This inhomogeneous, 

— iiV 

linear wave equation for ^h can be solved using the retarded scalar 
Green’s function for curved spacetime (DeWitt and Brehme 1960) : 


the retarded scalar Green's function for the curved \ 
( n „l)G(P T ,P) =| spacetime with the metric ( n -q)gp^ °f the (n,n-l) I 
approximation — a Green’s function with fractional J 
.errors £ n / 


(17) 


The result is 


i^ v (P) = 

n 


(n,n-l)[ l6,(- 8 )^ + ^ L )^ a , B **, a V (n-lf ^ 

( 18 ) 


4 , 
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This paragraph can be summarized as follows: Any approximation formalism 

of order (n,n-l), when augmented, by equation (18) for , , becomes an 
approximation formalism of order (n,n) . 

Special relativity and linearized theory provide a simple example 
of the above remarks: Special relativity is the approximation formalism of 
order (1,0) which one obtains by the extreme truncation process of setting 
h^ V =0 in equations (9) and (10): 




= 0 


= 0 

0 L-L 


O g pV ~ 71 pv * 


I 1 "''. 


V 


= 0 


(19 a) 


The retarded scalar Green’s function for a space with metric 
is 

0 G(P',P) = (4TT)" 1 « ret [f Tl pa (/ -x p, )(x CT -x a ')] . (19b) 

(Here <5 is zero if P lies in the causal past of P', and it is the 
ret 

Dirac delta function otherwise). Hence, equation (18) — by which one must 
augment special relativity in order to obtain a formalism of order (1,1) — 
has the form 


i^ v ■ / v < ? ') 5 re t'lv (xP - xP ' )( 


P - P W-x 5, )]dV 


= 4 


/ V (X°- I x-x’ I ,x’ ) 3 

— d x' 


( 20 ) 


x-x 


The resulting (1,1) formalism [equations (19) augmented by equation (20)] 
is the “linearized theory of gravity”; [see, e.g., §7.1, Box 7.1,. and. 
Chapter 18 of MTW.j 

Newtonian theory and the "quadrupole-moment formalism for wave gen- 
eration" are another example. Newtonian theory is the veak-field formalism 
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of order (2,1) which one obtains by not only truncating each series that 
appears in equations (9) and (10) , but by also imposing the slow-motion 
and small-stress assumptions 


2 _ 


0j 


v = |T J e^/T 


00,2 


T ij/ T °0| < 


( 21 ) 


(size of system)/ (characteristic time scale of changes) {.•.£ 


1/2 


Equation (18) , by which one augments Newtonian theory in order to obtain 
a formalism of order (2,2), has the form, when evaluated in the radiation 
zone 


— IT **TT 

2 h i j(t,x) = (2/r) I^Ct-r) = (gravitational radiation field) . (22) 

Here is the reduced quadrupole moment of the source, and TT 

denotes "transverse-traceless" part. This is the standard wave-generation 
formula of the quadrupole-moment formalism; see Chapter 36 of MTW. 

IV. WAVE GENERATION BY A WEAK-FIELD SYSTEM 
a) Motivation 

Weak-field systems are of two types: those with negligible self- 
gravitational forces (rotating laboratory rods; microwave cavities; etc.), 
and those whose internal motions are significantly influenced by self- 
gravity (pulsating stars; binary star systems; etc.). 

For a system with negligible self-gravity, special relativity gives 
a fairly accurate description of the internal motions; and, consequently, 
linearized theory [the (1,1) formalism obtained by attaching equation (18) 
or (20) onto special relativity] gives a fairly accurate description of 
gravitational-wave generation. 
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For most weak-field as trophysical systems, self-gravitational 

forces are important. In this case, when analyzing a system's internal 

motions, one must use a formalism of order (2,1); and when calculating 

.the waves those motions generate, one must augment the (2,1) formalism 

by equation (18), thereby raising its order to (2 ,2), 3 If the system has 

3 - * “ — — — 

in very special cases second-order gravitational forces may be as impor- 
tant, for the system's motions, as first-order forces. An example is a 
radially pulsating, weak-field star with adiabatic index very near 4/3 
(Chandrasekhar 1964); see also the discussion accompanying equations (61) 
below. When analyzing such systems one needs formalisms of order (3,2) 
and (3,3). 


slow internal motions and weak internal stresses, Newtonian theory [order 
(2,1) J will suffice for analyzing its motions, and the quadrupole-raoment 
formalism [order (2,2)] will suffice for wave generation. However, for 
analyzing fast-motion systems (e.g. , two stars flying past each other with 
hign velocity and deflecting each other slightly— the relativistic brems- 
st rah lung problem), one nee ds unrestricted (2,1) and (2,2) formalisms. The 
objective of the next two sections is to derive such formalisms. 

b) The Post-Linear Formalism 

A weak-field formalism of order (2,1), unrestricted by any constraints 
on velocities or stresses, can be obtained by truncating equations (9) and 
(10) at the appropriate order: 

^ - tT - 2 h^ 

(--jg) = 1 - L h = 1 + 1 'h 


(23a) 

(23b) 
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T av - r v T pa 

2 ,v l 1 av 2 1 1 av 2 X 


»“ 9 i*".„ - 2'“ V 


We shall refer to the formalism described by these equations as the 
"post-linear formalism". To analyze a system using the post-linear formal- 
ism, one must first specify the functional dependence of the stress-energy 
tensor 2 T yV on the system's nongravitational variables (e.g., density, 
pressure, velocities, electromagnetic field tensor, — ) and on the gravi- 
tational field 1 h UV ; and one must then solve equations (24a, b) simultane- 
ously for the system’s motions ( 2 T^ V accurate up to fractional errors e^) 

— uv 

and for the gravitational field (^h accurate up to fractional errors 
r \j £ ). Paper III will carry out such a calculation for the motion of two 
stars of arbitrary relative masses and velocities, which fly past each 
other with large impact parameter. 
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Having calculated a system's internal structure and motions using the 
post-linear formalism, one can then calculate the gravitational waves the 
system emits, 2 h^ V , evaluating expression (18). In evaluating (18) one 
needs an explicit expression for the retarded Green's function ;j G(P , ,P) as- 
sociated with the metric = Hp V + * . In the next subsection (§i) 

we derive ^G(P',P); then in iii we place constraints on our system which 
simplify ^G(P',P); and finally in §iii we use ^G(P',P) to evaluate the wave 
field 2 h PV . 



We shall obtain by taking the weak-field limit of the exact 

Green's function G(P',P) for a space described by an exact metric g 

• W 

The exact Green's function is formally rather simple, so long as the congru- 
ence of geodesics that emanate from the source point P* does not get 
focussed so strongly along the future light cone of. P’ that geodesics cross. 
Henceforth we shall assume "no crossing of geodesics on the light cone."' Later 
[eqs. (48), (48'), (48") below] we shall examine the constraints placed on 
the radiating system by this "no-crossing" assumption. 

DeWitt and Brehme (1960) have derived the exact Green's function 
G(P’,P) for the case of no crossing. Their Green's function consists of a 
"direct part" and a "tail" 


G(P',P) = G direct + G tail 


(25) 


The direct part is nonzero only if P lies on the future light cone of P f 
[denoted i + (P'>J. By virtue of the "no-crossing" assumption, when P is 

near j + (P') there is a unique geodesic from P* to P with a unique, squared 
length : • ; . 
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2 

ft f or timelike geodesic > rproper distance\ 

/•"World function 1 £or i m k Wlone geodesic ' 

n(P',P) = ( see SyI1 ge (1960) 1 2 U-I for spacelike geodesic alo g g 

(26) 

= o in notation of DeWitt and Brehme (1960). 

.+ 4 a w o = 0 G d;Lrect mus t have the form 

Because J (P') is characterized by S2 U , 

G direct (p , >p) = ( 4 lt ) _1 [A(F’,P)] 1/2 $ ret [m p ' > P)] 


(27) 


where 6 ret is the Dirac delta function on and near J + (P') , and is zero on 

and near the past light cone [i"(P')). ^he quantity A(P’,P) is an ampli- 
rude factor which would be unity in flat spacetime, but in curved spacetime 

is given by , 


A <P\P) 


dpt] | 3 Z «/3x a 3x B ' 1 1 
— — — 1,1 ’ , / 


|g(P) g( P ') 


1/2 


(28) 


f nr the tail different from, but equivalent 
We shall use an expression for the tan 

, _ _ Tn derive our expression we insert equa- 

te that given by DeWitt and Brehme. To deri 

tions (25) and (27) into the wave equation 


O G(P\P) = -[g(P)g(P , )]" 1/4 « x0 - x0 5 6(xl-xl 5 S(k2 ~ X } S(X } 

R V 


(29) 


The result is 


tail _ {( p A V2) 6(n)+[2 V Ki V 2 + (P s S!)A 1/2 ] 6' (A) 


□ G taix ■ -(4ir) 
s 


+ (Vf 2) 2 A 1/_ 6"(n)i , 


(30a) 


and second derivatives of the retarded Dirac 


where S' and S" are the first 
delta function, V is the 4-dimensional gradient operator, and 

variant derivative along the 4-vector 


is co- 
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K = W . (30b) 

[Here and below - we suppress Che subscript "ret" on b(Q).] We then manipu- 
late expression (30a) using the relations 

fl6 M (ft) = -26' (Q) , (vn) 2 = 2U , □ Q - 4 = -A -1 V A • (30c) 

[The first of these is a standard identity for Dirac delta functions; the 
second and third are eqs . (1.11) and (1.63) of DeWitC and Brehne (I960)]. 

The result is 

□ G tail = — (47T)™ 1 (0 A 1/2 ) S(Q) . (31). 

s s 

We then use relations (30c) and the relation A(P , P) = 1 to rewrite this 
in the form 

□ [G tail - (4ti) _1 (1-A 1/2 ) 6 ($2)] = +-C4tt) _1 07„ to A)i'(SJ) . (32) 

s * &- 

t^xl. 

Equation (31) tells us that G jumps from zero outside the light cone 
to a finite value inside the cone, without having any singularities on the 
cone. Equation (32) allows us to write (restoring the subscript "ret") 

3 if ?' £ I ~(P) (33a) 

(47f) _1 1 [to A(P’,p")] ja „ [fi(p',p")] ,a " 5; et mcp\p">] x 
X G(P",P)[-g(P")] 1/2 A" if P'SI(P) • (33b) 

Here I (P) means "the interior of the past light cone of P; T ' and condi- ■ 
tion (33a) suppresses the unwanted light-cone part of (33b) [i.e., sup- 
presses (47 t) ^ (1-A^ 2 ) 6 (ft)]. 
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Equation (33) is the form of the tail which we shall use. This form 

i 

was suggested to us by the work of Peters (1966). 

We now specialize the above equations for the retarded Green T s func- 
tion to the case of a weak gravitational field = ^liV + l^jJV ’ begin- 

ning with equation (26) for the world function. Let A be an affine param- 
eter along the geodesic linking P 1 to P 

ct 

C-(A) - geodesic with coordinates £ ( A) ; 

(3(0) = 3(1) « P|, ‘ 0$ A<1 . • (34) 

Then equation (26) can be rewritten in the form [cf. Synge (1960), p. 47] 

B(P',P) = I | g pv ( d ^ J/d A)( d ? V/ dA) dA . . . (35) 

3 

The right-hand side is actually an action principle for the geodesic equa- 
tion [cf. MTW, Box 13.3] . Therefore, if we evaluate the integral along the 
"straight line" 

0 3(A): C*(A) = + A(x a -x a ’) (36) 

(see Fig, 1) , which differs by a fractional amount of 0(e) from the true geo- 

2 

desic 3(A), we will make fractional errors in Q of 0(e ). Such errors are . 

2 

acceptable in ^G(P r ,P), since its frational errors are also 0(e ); cf., eq. 

(17). The result of integrating expression (35) along the slightly wrong curve 

q3( A) is 

i P-(P , ,P) = 0 W T ,P) + Y(P’,P) * (37) 

where 

0 Q(P’,P) = j X G X 2 * * * * * 8 n a3 , (38a) 
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Y (P\P) =| X a X 3 | ^dX , 
0 & 


v a cx a 

X = x - x 


(38b) 


(38c) 


Equation (37) is the desired expression for the world function. Turn 
next to the amplitude factor A(P'»P). Either by direct calculation from 
eqs . (28), (23b), (37), (38), and (B12), or by invoking eq. (95) on page 
63 of Synge (1960) , one arrives at the expression 


lA(P-,P) - ( 1 4i h ’li h,) det H^,aa’Il 

l 8 1 s 


(39) 


= 1 + 2a(P’ ,P) 


Here 


a (P',P) 5 i x a X 3 


(3 


l^a8 ^ (l - ^) » 


(40a) 


where is the Ricci tensor, accurate to first order in jh^: 


l R a6 


1 h f 

2 1 a8,p 


(40b) 


In eq. (39) we have simplified notation by using a prime to denote quanti- 
ties evaluated at P* , i.e., -jh' E ^(P 1 ) while ^h = jhCP) . Henceforth 
we shall reserve primes for this purpose— except that 5* and 6" are 
still derivatives of Dirac delta functions. 

a" 

Turn next to the "source term” (£n A ) , T ft* for the tail (eq. 33). 

>a 

The tail itself is of 0(e) compared to the direct part of the Green's 

function; therefore we can permit fractional errors of 0(e) in the tail — 

a" ’ 

which means we can use the zero-order value of Q* in the source of’ the 
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tail: 


- x a "-x a = X 01 ' 


(41) 


By combining this with equation (39) and by using equation (B7) of 
Appendix B, we bring the source of the tail into the form 


[in '^(P’.P")] „ [/!(?’, P")]’ a = B(P\P") , 


where 


6(P',P") ; j“T | ■ 

cA" 


dX 


(42) 


(43a) 


and is the "straight line" from P' to P" (see Fig. 1) 


: 9 a = x a ’ + Xx a " 


(43b) 


Turn, finally, to the propagator G(P",P) and the volume element 
1/2 A- 

(-g”) d x" which appear in equation (33) for the tail. Because the 
tail is of 0(e) compared to the direct part of the Green's function, we 
can ignore all curved space corrections in the amplitude of the propaga- 
tor (but not its phase) , and in the volume element: 


G(P",P)(-g") 1/2 d 4 x" = (4u) _1 S ret [ i n(P ,, > P)] dV 

in expression (33b) for ^G 


(44) 


tail 


All. of the pieces for the first-order Green's function are now at 
hand. By combining them [eqs. (25), (27) , (33) , (39), (42), and 
(44)] we obtain the following result: 


^(P* ,P)' = 1 G dlL ‘ eCt + 1 G tai1 


(45a) 
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(45b) 


^direct = (47r) -l [ i + a (P>,P)] 6 ret [ 1 S2CP',P)] . 

0 . if P' £ l"(P) 

(47T)" 2 | 8(P\P") «r e t [ i n ( p, ’ P ") ] 6 ret [ l fi(P " ,P)] ^ 

if P* c l“(P) * (45c) 

Here ^(P* ,P) [and similarly ,P") and ^(P",?)] is given by equa- 

tions (37) and (38), a(P‘,P) (the "focussing function”) is defined by ex- 
pressions (40), and &(?',?") (the "tail generator”) is defined by expres- 
sions (43) . 



Expression (45) for the Green's function is valid only if geodesics 
emanating from P T fail to cross on and near J + (P’). Crossing would be 
caused by gravitational focussing; and at any crossing point, the exact 
amplitude factor A(P' ,P) would diverge. Thus, the criterion for no cross- 
ing is finiteness of A along J + (P'). 

Consider our first-order expression (39) , (40) for^A. Evaluate it in 
the mean rest-frame of the source, with coordinates centered on P' so that 

X a - rn a , n°!= 1 , n = (unit spatial vector pointing from P T to P) , 
r = (spatial distance from P* to P) , 

X = r/r = (fractional distance from P 1 to P) ; (46) 

and invoke the first-order field equation = (j_ T a g " j_ T ) * Tfie 

result is 
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^(P'.P) = 1 + 2a 
r 

a = “ (n a n^ ^g) r(l-r/r.) dr . (47) 

O'. ‘ - . 

This expression for can never diverge if the source is bounded, because 

once the integration point r gets outside the source, ^T^g vanishes and 
^A stops increasing. However, if the focussing function ct approaches 
unity inside the source, then second-order and higher effects will come into 
play.. As one moves out into the vacuum beyond the source, those second-order 
effects will be essentially those of the "focussing" or "Raychaudhuri" 
equation; they will produce a divergence. Thus, the constraint 

CONSTRAINT: a(P’,P) « 1 for all P’ and P (48) 

is necessary for the validity of the first-order analysis, and simultane- 
ously protects us from "geodesic crossing." 

For a system that is roughly homogeneous with mass M and linear 
size L , equation (47) gives 

a ^ (m/l) ^ e <<■ 1 ; (48 *) 

so there is no problem in satisfying the constraint (48). However, for a 
highly inhomogeneous sytem (lumps of mass m and size & , separated by 
distances b » l ) , and for rays originating in one lump and passing through 
another, equation (47) gives 

a ^ b(m/£^H 'u (b/£) (m/£) . . 

.In this case the constraint' (48) is significant: it says that to avoid 
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too much ray focussing, the lumps must hot be too far apart 

(b/S,) « (JL/m) ^ lO 6 (£/R o )(M 0 /m) . (48”) 

.The Green’s function (45) would be much easier to use if, throughout 
it, we could replace the first-order world function ^£2 by its zero-order 
approximation = -j t \ X^X V . Let us examine [eqs. (37) and (38)] 

in the rest frame of our source, for points P on or near J + (P ! ): 

^(P'.P) - -|(x°+x)[-x°+ x + 2y(P',P)/(x 0 +x)3 

= X(-X°+X+ At ) (49a) 

D 

where 

X = ] xl = (distance from source to field point) (49b) 

Atg = y(P* ,P)/X = (’’Shapiro time delay") • (49c) 

For field points P far outside the source, the dominant contribution to 
the Shapiro time delay is the asymptotic "l/r" field of the source. It 
produces a huge delay of 

. _ /T v- r Shapiro time delay due to 1 

^asymptotic field of sourcej 

M = (mass of source) (50) 

L E (characteristic size of source) 

This delay is time independent and is independent of where inside the 
source P T is located (aside from a negligible piece of size 'v 2ML/X) ; 
therefore its only effect on the radiation is to delay the arrival time 
at a given radius. Henceforth, for ease of calculation, we shall remove 
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this constant delay from the argument of our Green’s function. We can 
always reinsert it at the end of the calculation if we wish. With this 
constant delay removed, we can rewrite aS 

CUP 1 ,P) « X[~X° + X + (Atg- A)p T p] 

= 0 n(P*,P) + Cy (P’.P) + ax\], (5] 

for P’ inside the source and P far outside it , 



nj m £n(b/£) for lumpy source 

Henceforth we shall assume that this internal time delay is small com- 
pared to the characteristic timescale on which the source changes i.e. , 
small compared to the characteristic reduced wavelength * of. the radia^ 

tion emitted, 

CONSTRAINT: (At s -) interna i ^ m £n(b/£) « * * 

[Example: If -k is 100 times larger than the Schwarzs chi Id radius, 2m, 

4 

of a lump, then b/£ can be as large as exp (10) ^ 2 x 10 without 
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causing problems. Another example: If b (which is the case for brems- 

strahlung) , and if. & » m (which is required for fields to be weak), then 
the condition b ^ £ (separation of lumps bigger than size of lumps) guarantees 
that constraint (54) is satisfied.] 

The constraint (54) allows us to expand our delta functions 6 (^) 

in powers of the internal time delay. Discarding terms that are quadratic. 

and higher-order in (At_) . _7^, we obtain for the Green's function (45) 

*> internal 

.jGCP'.P) = iG dlreCt + 1 G tail , (55a) 

/ 

1 G dlrect (P',P) = {6 ret (-j X a x 6 n a g) + a(P*,P) 8 ret (| rf n a g) 

+ [y(P',P) + AX W U u ] 8; et (ix a X 5 n a3 )} . ; (55b) 


tail(p 1)P ) = , 


0 if ?' X"(P) 

-2 


-w 


e(P’,P") V 3 ret (| 3 ? t 5 e n„ R )d 4 x " 


aB' 


if P € T(P) . 


In these equations 


v a _ a a 

X = x - x 


_ a .. a" 

X - x - x 


„a" _ a" a 1 
X = x - x 


(56) 


see figure 1, 

Equations (55) are our final form for the scalar Green's function 
in a space with linearized metric -jg ~ H + l^uv * Green's func- 


tion has fractional errors 


| (errors in ^ Maximum of {e ,oe , [ ( At s ) ittternal /^] £ } iri general, 

2 

'u e for most sources ; ( 57 ) 

and it has been stripped of its asymptotic time delay (eq. 50). 
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ill. The Gravitational-Wave Field 2 “ 


By inserting expressions (55) for ^G(P’ ,P) into equation (18) we ob- 
tain the following expression for the gravitational field far outside a 
weak-field source: 


2^W ' 2“TL 




h^ V = 4 
2D 


ret'2“ M a3 2 / L 1 

,h^ V = 4 | a (P\P) n aB ) 2 T PV (P’) dV , 


h^ v = 4 
2 n TR 


V retv 


! ctS ; 2 



(58a) 

")] dV , 

(58b) 

.V 

(58c) 

T^CP') d 4 x’ , 

(58d) 


-niv 

2^ 


2 n TL 


4 j S ret (~^ a3 > L x Cl + l^a 

(— 1/tf) 


A , 

a x 1 

at ?’ (58e) 


SCP’.P") «; et (f^V\ e ) S ret (|r¥ n„ fl ) 


aB' 


P’ €l-(P) 


X 2 T UV (-P t ) d 4 x M dV 


(58f) 


a T a" 

Here P' and P" are source points with coordinates x and x (cr.fig.l); t 

field point P has coordinates x 01 ; Q T^ V , - t^ V T , and -h.^ are the stress- 
energy tensor, the pseudotensor, and the gravitational field obtained by a 
post-linear analysis [eqs. (23) and (24)]; 6 is the Dirac delta function 01 

(X 

the future light cone of the source and zero on the past light cone; X , 

■pt 


TpJ. 1 ir 01 

X , and X are 


'a _ a a 1 pi _ a a" „a" _ a” 1 

X — x — x , X = x - x , X -x -x 


.a” _ a” a r 

— V - X 


(59a) 


a, 3, and y are defined by integrals of the first-order Ricci tensor 

and of the metric perturbation ^h^ v along the straight line between 


two points 
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a (P\P) = “X a x^ 


^^(x^ + AX P ) A(1~A) dX 


rt” ft*' 

3(P’,P") = x a x p 


n R D (x y, + XX U ") X 2 dX 

1 dip 


•1 


Y(P’,P)=I x a X £ 


f AX^) dA 


(59b) 

(59c) 

(59d) 


-AX a U^ is that portion of y which is produced by the asymptotic, 1/r, 
external field of the source 


. ct __ . ,_ x a. \ / Shapiro time delay \ y { distance from source^ 

a a \ produced outside source/ \ point to field point/ . 

(59e) 

(see §ii above); and P T € I (P) means that the integration (58f) is per- 
formed over field points P’ that lie inside but not on the past light cone 

of P . 

Each piece of the distant gravitational field 2 h“ V has its own 
physical origin and significance; 

2h^ V is the "direct field." It is produced by the stress-energy 
T yV and propagates as though spacetime were flat. It includes the zero- 
order, non-radiative , "l/r" field of the source, and also that portion, of 
the radiation produced "directly" by the source's motions. If the internal 
gravity of the source has negligible influence on the source's structure and 
evolution, then all other parts of 2 ^^ negligible compared to the 

direct field ("linearized theory"; cf. eq, (20) and the associated discus- 
sion) . 

2h yV is the "focussing field". It is the amount by which the direct 
F 

field is augmented due to focussing as it passes through regions of nonzero 


Ricci curvature (nonzero stress-energy) . 
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_h^ is the ’’transition field” [first discovered in the equations 
Z TR 

of general relativity by Chitre, Price, and Sandberg (1973,1975); analogue 
of "electromagnetic transition radiation", Ginzburg and Frank (1946)] , 

It is the amount by which the direct field changes due to Shapiro-type 
time delays within the time-varying source. 


2 *^ is the "whump field." It is the field generated by "gravi- 
UV — 1 "ijo 

tational stresses” ^t£ ^ + (16 tt) ^h a 1^ p ' k ave given it the 
name "whump" because in our minds we have a heuristic image of gravita- 
tional stresses linking various pieces of the source, and going "whumpity- 
whump-whump" as the source contorts and gyrates. 


is the "tail field",- It is generated by the direct field in 
those regions where focussing has deformed the geometry of the direct 
wave fronts. 


Although it is useful, heuristically and in calculations, to split 
2 ^^ into these five pieces, one should not attribute too much physical 
significance to each individual piece. For example, no individual piece 
satisfies the Einstein field equations or the deDonder. gauge condition. 
However, the five individual pieces combine in such a way that their sum 
does satisfy the field equations and gauge condition; see Appendix C. 
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V. S LOW-MOT ION LIMIT OF THE WAVE -GENERATION FORMULAS 

Consider a weak-field system which has slow internal motions and 
weak internal stresses. Characterize it by the following parameters: 


L = (size of system) 


* = (characteristic time-scale of system) - 


M H (mass of system) 


(reduced wavelength 
of radiation) 


(60) 


v 


- (It^I/T^) = (maximum internal velocity) 

v 1 1 ' max 


S 2 = (| T 1 ^ 1 /T 00 ) - maximum of (stress) / (density) 


Chapter 36 of MTW derives the quadrupole-moment formalism for . gravitation- 
al wave generation under the following assumptions (eqs. 36.18 or MTW) . 


L/)c « 1 which implies v « 1 i 


M/L « h/k , S « L/A 


(61a) 

(61b) 


Constraint (61a) is the standard slow-motion assumption— the only assump- 
tion truly necessary for validity of the quadrupole-moment formalism (see 
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Paper IX ). Constraints (61b) say that the motion must not be too slow' 
i£ a we ah- field calculation is to yield the quadrupole-moment formalism. 
In terms of the characteristic frequency “ = 1/ *■ this not too 
assumption says 


to 2 » (M/L)(M/L 3 ) , u 2 » S 2 (S/L ) 2 


( 61 b') 


A violation of these assumptions occurs, in dynamical systems, only when 

the gravitational and stress forces counterbalance each -other so precisely 

that second-order gravity, can affect the motion significantly [ cf. 

Chandrasekhar (1964)]. In this case an analysis based on the post-linear 

1 - oive a correct description of the radia- 

approximation cannot possibly give 

tion. 

It is instructive to see how the post-linear radiation formulas (58) 
of this paper yield the quadrupole-moment formalism, when applied to a 

system satisfying constraints (61), 

We begin by combining the direct and whump fields (58b, e) and ' then 
breaking them up again, differently , 


t-uv, + h^ v ; 

2 h D + 2^ ~ 2 b DWl 2 h CW2 


h^, 5 4 
2 DW1 


(62 a) 


(62b) 


^et^ 0 V [C "1 S) ( 2 TVV+ l t L-L )] a t P’ d * x ' ’ 

, AviyS .r rUP ] dV • (62c) 

5 ret ( 2 X ^ ^asll ,a l n ,P J at P' 


& 2 


-rest frame of the source 

We then evaluate expression (62b) in the rest tram 

(- 1 g)(? T ^ ; + , l t L-I?* e t d 3 x t ; 


= 4 

2 mi 


(63) 


|x - x 
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and by carrying out the analysis of^'MTW §36.10, we bring the spatial 

State that (-^g) ( ? T ;iV + here plays the same role as T uv + t^ v 

in MTW. The key properties which they share are (i) vanishing coordinate diver- 
gence; (ii) same role in retarded integral for ^h 


transverse-traceless part of this field into the form 

[^ 1 (t,x)] TT = (2/r) (d 2 /dt 2 ) i^(t-r) 

2 (64) 

^ (M/r) (L/*) z 

TT . 

Here is the "reduced quadrupole moment" of the source, and 1S lts 
transverse traceless part. This is the standard quadrup ole-moment formula 
for the radiation field. • ■ 

An order-of-magnitude analysis shows that all other parts of our 

expression (58) for 9 h , are negligible. In particular, by using the fol- 

lowing relations valid for the source's interior 


1 h°° ^ m/L , .jh 03 * ^ Mv/L , 'v MS 2 /L, T 3 ’ k = MS 2 /L 3 , 

^ Q * nj ^ /h 9 3 * M/L ’ (65a) 


as well as the relations 


a ^ M/L , (y + ATiy % rM , 


(65b) 


we obtain for the ratio of each other part to the "DWl" part (eq. 64): 


r-jk ,, r~jk . TT i 11, f S ,2 . 

2 h DW2 /( 2 h DWl ) 'I " + L/-^ 


(66a) 


rjk, , Hk v TT | 
2 h F / 2 h DWl 1 




(66b) 
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(66c) 


7“j k , . r*j k N TT 
2^TR^2 h DWl^ 


0/ 


P k /( h^ k 

2 TL ^2 DW1 


TT i M S 

) I ^ (r)( 


T) ( l7X } k< 1 


VI. SUMMARY OF OUR "PLUG-IN-AND-GRIND" FORMALISM FOR WAVE GENERATION 

Our post-linear formalism for wave generation can be summarized as 
follows: 

Regime of Validity . The formalism is valid for any system satisfy- 
ing these constraints: (i) The gravitational field must be weak everywhere 

. _jv vl I 

|h H I « 1 everywhere , (67a) 

and the source must be isolated [see discussion preceding eq. (4)].. 

(ii) Gravitational and nongravitational forces must not balance each other 
so precisely as to enable second-order gravity to influence the system's 
motions significantly, (iii) The source must not focus substantially light 
rays emitted from within itself. Mathematically this constraint says 

|a(P',P)| « 1 for P' any event inside the source, 

P any event on the future light cone of P T , 

(67b) 

where a is defined by equation (59b). For further discussion of this 
constraint, see the first half of §IV.C.ii. (iv) The "Shapiro time delay" 
for light propagation within the source must be small compared to the 
characteristic timescale ^ for internal motions of the source. Mathemati- 
cally this constraint says that in the mean rest-frame of the source 

(At S> internal = y(F’ ,P) / 1 x-x' | « * . (67c) 

Here x and x 1 are spatial locations of events P and P' that lie 
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inside the source, P is on the future light cone of P r , and Y is de- 
finedby equation (5 9d) . For further .discussion, see the second half of 
§IV .C.ii. 

Calculation of the System's Motion . For a system satisfying these 
constraints one calculates the internal structure and dynamics by using 
the post-linear formalism of §IV.B[eqs. (23) and (24)]. 

Calculation of the Distant Field . To calculate the gravitational 
field 2 h^ V in the radiation zone, far from the source, one takes the re- 
sult of the post-linear analysis, plugs it into equations (58) and (59), 
and grinds . 

In paper III we shall use this formalism to calculate gravitational 
bremsstrahlung radiation. 
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APPENDIX' A 


■ WHY USE THE CURVED-SPACE SCALAR-WAVE OPERATOR? 

In laying the foundations of our analysis [in and near eq. (6)] we 
write the Einstein field equations in terms of the curved-space scalar wave 
operator 0 • We choose to do this because the obvious alternatives (the 
flat-space wave operator or the curved-space tensor wave operator D^.) 
would ultimately lead to complications or dangers in our analysis. 

The flat-space operator D f treats the field propagation from the out- 
set as though it were on flat-space characteristics (straight coordinate 
lines) . Because the true characteristics suffer the Shapiro time delay 

which involves a logarithm of distance, the use of would lead to logar- 
in the radiative field at large r . 
ithmic divergences^ If one were sufficiently careful, one could remove 

those divergences without serious error — but that is a dangerous enterprise. 

Even if one succeeded, one would be left in the end with the interesting 

effects of focussing, time delay ("transition radiation") , and tail all 

lumped into the "whump" part of the field. We prefer to keep them separate. 

Consider next the curved-space tensor wave operator 

□ P V 5 i^ V ;a + 2R V - 2R (A1> 

t a a 3 a 

[cf. MTW eq. (35.64)]. Because the true propagation equation for very weak 
gravitational waves on a curved background is = 0 , it is tempting to 

formulate our analysis in terns of □. rather than □ . . By using □ we push . 

t S 5 

into the "whump" part of 2 h ^ V an important physical effect: the curvature- 

induced rotation of polarization. In effect, part of our whunrp field corrects 

the error in our direct field's unrotated polarization. Had we used 

rather than □ , polarization rotation would have shown up in §IV.C.iii as a 
s 


37 


separate piece of the radiation field. 

The tensor wave operator has a disadvantage which, for our purposes, 
outweighs the above advantage. Suppose that one constructed a tensor 
Green’s function for 


□ G^ a ' S ’(P',P) - 


f(g- 


g VS ' + 8^' 


va\ , , s -1/4 , , ts 

g )(gg ). 6 . (x-x ) 


CA2) 


or for any other wave operator with the form 

h^ V = h^ V;a + (any "background" field)^ V h'^ . (A3) 

other a 

That Green’s function would have a first-order tail 1 G Cail with 

"sources" ^ involving the Riemann tensor [cf. eqs . (43a) and (45c)]. 

Such a tail would originate everywhere on the light cone of P* , whereas 
the tail 1 G tai1 for our scalar Green’s function originates only on rays 
that have passed through matter. In practical calculations involving 
lumpy sources — ’see, e.g. , Paper. Ill — that tail would be as difficult to 
calculate as the whump part of the field. We prefer our scalar tail be- 
cause of its greater simplicity. By using □ we dump all serious calcu- 
lational complexities, for lumpy sources, into the whump part of the field 
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APPENDIX B 


LINE- INTEGRAL IDENTITIES 

The weak- fie Id Green’s function -^(P^P) used in this paper involves 
three integrals a v S,y along "straight lines". In this appendix we take 
the line of integration to be 


H x“' + AX 3 , 0 < X < 1 , x“ = x“ - x“ 


(El) 


The three line integrals are 

Y | xV | x h uv -dA 

0 
• 1 


a = | xV f X(X-X) .dX 


B = xV 


;2 


Av x dX ’ 


-o 


where ^ is assumed to satisfy the deDonder condition 


(B2) 


(B3) 


(B4) 


l h ViV ,V “ l h UV ,V “ 2 l\u " 0 


(B5) 


and the Ricci tensor is therefore given by 


* - _ I h p 

l K U v 2 1 UV.p J 


(B6) 


and where the index notation used is that of a Lorentz frame in flat spacetime 
Below we list a number of useful identities linking the line inte- 
grals Ci»3,Y. , their derivatives at point P , and the values of 

and . R at P • 

1 yv 


xP “,p = i s 


(B7) 
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d a 
X K X a 


,pa 


_ 6 + I xV l R pa 


(B8) 


X P Y • = y + ~ x p x a h 

' ,p 1 2 1 pa 


(E9) 


x p x a Y ipa = x p x a lhpo + j XW lhpa>T , 


(BIO) 


y 


p = -6 + .h 


(Bll) 


Similar identities involving derivatives at P r and mixed derivatives at 
P* and P can be derived fairly easily. For example. 


Y 


P 


P’ = -2a -|l h -il h ' 


. (B12) 


where = ^h(P) and ^h’ = ^h(P*) . • . . 

The derivations of these identities are quite straightforward. The 

necessary techniques are illustrated by the following derivation of iden- 
tity (B7): By differentiating definition (B3) and making use of equations 
(Bl) , we obtain 

1 ' ■ 


a * x 

*P 


V 


R A(l-X) dA 
PP 


0 


+ 1 x y x v 


(3R^ v 7aC a )(3S a /3x P ) A(l-A) dA 


= x y 


R A(l-A) dA + ~ X U X V 

yp 2 


- X y [ R A(l-A) dA + ^ X y X V 

1 yp 2 


( w (X6 V xa “ x) dx 


R A (1-A) dA 
yv,p 


Here R 


yv,a 


8r^^/ 9^ G is the derivative of R^ at the integration point 


(3(A) . When multiplied by X p this expression gives 
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X p a = X y X V 
»P 


R pv A(l-A) dA +~xV j (R yv>p X p ) A 2 (l-A) dA 

' o 
1 


- X y X V j A(l-A) dA + ~ X y X V j (dR^/dA) A 2 (l-A) dA . 


0 


0 


By integrating the last expression by parts we obtain 

1 

X p a = X U X V f RJX(I-X) - ~(d/ d A)(A 2 - A 3 )] dA 
»P ) ^ 


0 


- 1 *v 


V * 2 dA = i e 


QSD. 


In this case the integration by parts gave no endpoint terms; but In other 
cases [eqs. (B8)-(B12)] nonzero endpoint terms are obtained. 

In manipulations of our weak-field Green’s function ^G(.P’,P) and. of 
our second-order gravitational field ,-,h^ V (see, e.g. s Appendix C) two other 
identities are useful: 


86’ = (a6) p - a p 6 . + 4n ao.(x-x') , 

>P »P ^ 


(B13) 


(v<5’+ ad) p =a P 6 + 6 

»P »P 1 » 




pa 


-[a + | x h + (Y^ D X Q 1 h ao )3 P J 4tt 6 4 (x-x’) 


I x a 

,p 2 A rap' 


(B14) 


Here 5 is the flat-space propagator between P ’ and P 


656 fi X P X CT n ) 
ret 2 • pa 


(B15a) 


which is related to the 4~dimensional Dirac delta function 
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(B15b) 


<S 4 (x-x t ) 5 6 (x°- x°') Six 1 - x 1 r ) 6(x 2 - x 2 ') d(x 3 - x 3 ’) 

b y 

(S ,p P = ~ 47T ; (Bi5c) 

and $ T is the derivative of the propagator (B15a) with respect to its argument 
The absence of primes on indices and onh's in (B 13) and (B 14) indicates that 

all derivatives and endpoint terms are taken at P ; none are at P* . The 
identities (B13) and (B14) can be derived with some labor from the iden- 
tities (B7)-(B11). 
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APPENDIX C 



PROOF THAT THE "PLUG- IN-AND- GRIND " FORMULAS FOR SATISFY' 

THE FIELD EQUATIONS AND GAUGE CONDITION 

Here we briefly sketch the proof that our second-order gravitational 

field [eqs. (58)] satisfies the second-order Einstein field equation [eqs. (16) 

with n = 2] and the deDonder gauge condition ? h^ V =0 . As part of our 

z . ,v 

proof we shall derive expressions for the amount by which each piece of 
^h^ fails, by itself, to satisfy the field equations and gauge condition. 

A preliminary step in our proof is to rewrite the "tail" and "transi- 
tion" fields (58f) and (58d) in new forms. 

Although expression (58f) for the tail seems optimal for practical 
radiation calculations, the restriction P’ e I (P) makes it nasty for 
formal manipulations. To get rid of this restriction, we take expression 
(B13) for 36' , in it we replace p by P'\ and. then we insert it into 
expression (58f ) . The result. 


2^L = (1/ *> ft V Wl ^ W 


* 2 X^ V (P * ) d 4 x" d 4 x r 


is an expression which gives the same value for 2 ^TL w ^ iet ^ er one imposes 
or omits the restriction P T e I (?) . One way to see that (Cl) is oblxvi- . 
ous to the restriction p’ € I~(P) is this: Take the source equation (31) 
for the tail of the exact curved-space Green's function; calculate its lowest- 
order form 

G tail j0 P = -(4TT)" 1 [a(P t ,P)] jp p 6 (~ x a x 3 n ajS ) ; 
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invert this using a flat-space propagator; use the resulting G to 
calculate 2 h^ ; the result will be expression (Cl)— and nowhere in the 
derivation did one need to impose the restriction P r € I (P). 

Expression (58d) for the transition field involves a "time-delay 
function" yi?' ,P) from which the logarithmic, "external time delay" 

A(-X a U a ) has been removed. A straightforward subtraction of the external 
time delay is well suited to practical calculations, but poorly suited to 
formal manipulations of 2 ^^ * In formal manipulations of this appen- 
dix we shall perform the truncation in a "smoother" manner: We surround 
the source by a (hypothetical) cloud of negative-mass material, with total 
mass, -M, equal in magnitude to that of the source, -HI. We put the cloud 
far enough from the source (e.g. , at radius s t % 100L) that it is very dif- 
fuse, and this contributes negligibly to the line integrals a and 3 ; . 
but near enough that the Shapiro time delay 2M £n(£/L) in going from 
source L to cloud £ is small compared to the timescale A of the 
source’s internal motions. The cloud automatically removes the external c; 
Shapiro time delay; no artificial truncation of y is needed. The second- 
order gravitational field is then given by equations (58) and (Cl) every- 
where (inside the source and out), except that we must remove the artificial 
truncation from (58d): 

= 4 | y(P',P) 5; et (f X a X® n a0 ) 2 T UV ( P ') d 4 x' . (C2) 

Turn now to the proof that our second-order field satisfies the 
second-order Einstein field equation. We begin by applying the first-order 
wave operator 

1 n s = (n a3 - h a3 ) . • (C3> 
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to each of the five pieces of our second-order field. 3y applying 
to the direct field (eq. (58b)) and by using equation (315c) we obtain 


ft 2^D 16lr ^ 2 l 1 ^ 2 T ^ 


4h^ a 3 3 . 

P o- J 


5 ret ( i rf -W ^ 


(C4a) 


By applying , □ to the whump field (eq. 58e) and by using (315c) we 

x S 

obtain 


n ryv rVP u va 

ft 2 w = " l6Tr 1 t L-L ! h ,a l h .,p 


(C4b) 


By applying 1 0 to the tail field ( eq. Cl) and by using (B15c) we obtain 
J- s 


1°. 2^L - - 4 [“0”.™/ 5 ret^ xV V dV * (C4c) 


v CD » \ r 


By applying □ to the focussing field (b8c) , and by using (B15c) and the 

-L S 

relation a(P,P) = 0 [cf . eq. (B3)] we obtain - 


l°s 2<= 4 f ta(P’.P)] >p p 6 ret (|x“x s n a g) / V CP’) dV 


(C4d) 


+ 8 


[ct (P 1 ,p>] >P [S ret c| x a x 0 n a g)] ,p 2 t’ jv (P’ ) d 4 x’ . 


By applying to the transition field (02), and by using (314), (Bll) , 

J*. s 

(B15c) , and limits as P' -> P that are obtainable from (Bl)-(B4),we obtain 


1° 


5 ret<t X “ xB V 2 tUV < P '> 


(C4e) 


- 8 


[a(P’,P)] „ [S.„(Vx¥ n„)]’ P ,T !JV (P') dV , 


,p ret 2 


ct0‘ 
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By adding up all five pieces (C4a)-(C4e) we obtain 

1-t Ti-iV 1 Tn m'dV , UV -i rllp TyVa 

l D s 2 h ” 16 '‘ I(1 2 l h) 2 T + l t L~L l h l h ,| 


(C5) 


which is the second-order Einstein field equation (16) . 

Turn now to a proof that our field (58) satisfies the deDonder gauge 

condition 2 ^ v = ^ except for fractional errors of O(e^). From (58b) and 

the relation 
r 
3 


v 


6 ret (| x“x e n aS ) f (P’> A' - 


^ Stijft 6 ti ae ) 3 v ,f(P') dV (C6) 


valid for any function f(P'), we obtain 


2 h D V ,v“ S ret ( I tl- ^CP')]}^, dV . (C7a) 


From (58e) , (C6) » and a = 0 we °^ ta ^ n 


2^ V a V = 4 


‘Wl 5 ^ V + l^.av l^.phtP- d 


4 x’ 


(C7b) 

We now add (C7a) and (C7b) and use the post-linear equations of motion (24a) 
rewritten in the form 




[cf. eq. (3)] to obtain 


<2^ + 2^\v " V l^'Vv’ iT <?"> ,p- dV 


P ' 

We then use an integration by parts on x together with (C6) and a re- 
labelling of indices to obtain 
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, Apv 

2 n D 


+ 


2\ ) 


A 


= (4tt) 


‘ (5 (-L r» N h^^CPM ^(V") 

°ret ( 2 XX n aB ) l h V'p' l h - C } 


v ‘ 


dV 


We then give P’ the new name p” and rewrite h^CP") as a retarded 
integral [the solution to eq.. (24b) }; the result is 


<2*T + 2<>,v ■ (1/W) 3 V \\ 5 *,A^ W l^ a(P "> .* 


x [6 ret (i X a "x 3 \ e )] )p „ a „ 2 A<P')dV A” . >7c) 


By applying 3^ to expression (Cl), adding it onto (C7c), using identi- 
ties (B14) and (B15C), integrating by parts, and using the limiting forms of 
(B1)-(B4) as . V -*■ P r , we' obtain 


/ AiV , 
( 2 h D 





h^ V ) 

2 TL J , v 


Y(P’,P) 


a(P',P) « tet (fxV n ae ) 


2 T VV (P') d 4 x' 
2 # V (P') d 4 x' 


(C?d) 


Comparison with expressions (C2) for 2 h)^ and (58c) for 2 h^ V shows 
that 


VC + 2^" 


h^ V 
2 TL 


■ruv , tij\x 
2 n TR 2 h F ) , v 


= 0 


(C8) 


i.e., our total second-order field does satisfy the deDcmder gauge con- 
dition. 
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FIGURE CAPTIONS 


1» The points ? , P" used in evaluating the post-linear 

Green’s function ^G(P* ,P) and in calculating the post-linear 
gravitational-wave field ^h^CP) . Part (a) shows the parametrized 
straight-line curves Q C(A) q^ 5*(A) linking P, P r } and P"; 

11 ] J ^ 1 J 

part (b) shows the 4 -vectors X , X , and linking them. 
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